Gravitational constant calculation methodologies 
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We consider the gravitational constant calculation methodologies for a rectangular block of the 
torsion balance body presented in the papers Phys. Rev. Lett. 102, 240801 (2009) and Phys.Rev. 
D. 82, 022001 (2010). We have established the influence of non-equilibrium gas flows on the obtained 
values of G. 
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The methodologies of calculating G are described in 
[1]. In methodology 1, the periods of anharmonic oscil- 
lations are determined by the Runge-Kutta method. In 
methodology 2, calculations are carried out by analytic 
formulas after expanding the attraction torques in series 
in odd powers of the balance deflection angles f. Compli- 
cated shapes of the bodies prevent obtaining an analytic 
representation of the attraction torques. In [2-4], the at- 
tracting bodies have a spherical shape while the working 
body is fabricated as a rectangular quartz block coated 
with two thin metal layers. In the model system, the 
block is changed for a thin rod of the same mass. In- 
creasing the distance between the rotation axis and the 
attracting bodies provided an adequacy of the model and 
real systems at small if. A three-position scheme of mea- 
surements has been considered. In position 1, the attract- 
ing bodies are placed at the ends of the block, in position 
2 they are absent, and in position 3 they are rotated by 
7r/2. The positions, oscillation periods and values of Gij 
according to methodologies 2 [1] and 3 [2-4] for one in- 
verse measurement run (positions 3, 2, 1) and one direct 
run (positions 1, 2, 3) in the first experiment as well as 
one direct and one inverse run in the second experiment 
are given in the table. The results are actually indistin- 
guishable. Methodology 3 understates the values of G in 
combinations 1-3 and 3-1 by 6 ppm, while in combina- 
tions 2-3 and 3-2 it slightly overstates them by a small 
nonlinearity persistence in the equations of motion. A 
correction of 212 ppm for inelasticity of the suspension 
thread was not taken into account. It contradicts the 
dislocation theory of inner friction [5] . After removing a 
calculation error made in [1] and using the actual values 
attraction torgue data from [4], average value of G as 
compared with [2] in first experiment was overstated by 
192 ppm and in second experiment was understated by 
206 ppm. In this case the first five values of Gij under- 
state the average value of G at 1120 ppm, and overstate 
its last eight to 386 ppm. Only this eight have to be taken 
into account, since they reduced the standard deviation 
of 24 times. The difference in values of G for direct and 
inverse runs was in the first experiment 1016 ppm, and 


the second one 472 ppm. A monotone drift of the oscil- 
lation period, the differences in Gij in direct and inverse 
runs, appreciable deflections of Gij from the normal val- 
ues when taking into account position 2 indicate the exis- 
tence of slowly decaying non-equilibrium gas flows in the 
chamber. An unlucky choice of the material and shape of 
the working body has strengthened their effect. It could 
be weakened by increasing the density of the block and 
the attracting masses. The equality of all combinations 
of Gij was provided by an oscillation period diminished 
by 15 ms in the flrst experiment and increased by nearly 
5 ms in the second experiment. One also cannot exclude 
the influence of a magnetic interaction [1], which is hard 
to single out against a more powerful factor. 
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